Exact exponents for the spin quantum Hall transition 
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We consider the spin quantum Hall transition which may occur in disordered superconductors with 
unbroken SU(2) spin-rotation symmetry but broken time-reversal symmetry. Using supersymmetry, 
we map a model for this transition onto the two-dimensional percolation problem. The anisotropic 
limit is an sl(2|l) supersymmetric spin chain. The mapping gives exact values for critical exponents 
associated with disorder-averages of several observables in good agreement with recent numerical 
results. 

PACS numbers: 73.40.Hm, 73.20.Fz, 72.15.Rn 

Noninteracting electrons with disorder, and the en- 
suing metal-insulator transitions, have been studied for 
several decades, and are usually divided into just three 
classes by symmetry considerations. Recently, the ideas 
have been extended to quasiparticles in disordered su- 
perconductors, for which the particle number is not con- 
served at the mean field level. Several more symmetry 
classes have been found [[[]. One of these, denoted class C 
in Rcf . jy , is of particular interest || tl ■ This is the case 
in which time-reversal symmetry is broken but global 
SU(2) spin-rotation symmetry is not, and spin transport 
can be studied. In two dimensions (2D) it can occur in 
d-wave superconductors. Within class C, a derealization 
transition is possible in which the quantized Hall conduc- 
tivity for spin changes by two units, resembling the usual 
quantum Hall (QH) transition but in a different univer- 
sality class. When a Zeeman term is introduced which 
breaks the SU(2) symmetry down to U(l), the transition 
splits into two that are each in the usual QH universality 
class. 

In this paper we present exact results for a recent model 
Mm for the spin QH transition, in class C, in a system of 
noninteracting quasiparticles in 2D. We use a supersym- 
metry (SUSY) representation of such models, considered 
previously |q] , to obtain a mapping onto the 2D classical 
bond percolation transition, from which we obtain three 
independent critical exponents, and universal ratios, ex- 
actly. An anisotropic version of the model is also mapped 
onto an antiferromagnetic sl(2|l) SUSY |7]] quantum spin 
chain. The results are in very good agreement with re- 
cent numerical simulations |Jg . 

We study the spin QH transition in an alternative de- 
scription that is obtained from the superconductor after 
a particle-hole transformation on the down-spin parti- 
cles |2j, which interchanges the roles of particle num- 
ber and z-component of spin, and so particle number 
is conserved rather than spin. This makes it possible to 
use a single-particle description, at the cost of obscur- 
ing the SU(2) symmetry. The single-particle energy (E) 
spectrum has a particle-hole symmetry |jj under which 
E — > —E, so when states are filled up to E — 0, the 



FIG. 1. The network model. 



positive-energy particle and hole excitations become dou- 
blets of the global SU(2) symmetry. In this picture, a 
(nonrandom) Zeeman term H for the quasiparticles maps 
onto a simple shift in the Fermi energy to E = H |pj, 
splitting the degeneracy. 

The model [f||| is a network (generalizing Ref. 0), in 
which a particle of either spin and with E = 0, repre- 
sented by a doublet of complex fluxes, can propagate in 
one direction along each link (Fig. pi). The propagation 
on each link is described by a random SU(2) scattering 
matrix (the black dot), with a uniform distribution over 
the SU(2) group; the absence of an additional random 
U(l) phase here is crucial and implies that the global 
SU(2) spin-rotation (or particle-hole) symmetry is unbro- 
ken. As in Ref. M, there are two sublattices, A and B, 
on which the nodes are related by a 90° rotation. Scat- 
tering of the fluxes at the nodes (black squares) is de- 
scribed by orthogonal matrices diagonal in spin indices: 
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where S — A, B labels the sublattice and a =|, J, la- 
bels the spin direction. The network is spatially isotropic 
when the scattering amplitudes on the two sublattices are 
related by t 2 Aa + t| CT = 1. 

The network has a multicritical point at tAa = tsa = 
2 -1 / 2 (for the isotropic case). Taking tAa 7^ ^Bo (but 
keeping ts-\ — isj.) drives the system through a QH tran- 
sition between an insulator and a QH state, and the Hall 



conductance (now for charge) jumps from zero to 2. Mak- 
ing ts-\ ^ tsi breaks the global SU(2) symmetry, and 
splits the transition into two ordinary QH transitions Q 
each in the unitary class. As we will argue later, this 
perturbation is different from the uniform Zeeman term. 

We briefly describe, for the present case, the main steps 
of the SUSY method for the network models M . Trans- 
port and other properties of the network, such as its 
conductance, may be expressed in terms of sums over 
paths on the network. Such a sum may be written 
in second-quantized language as a correlation function, 
(...) = STr (T. . . U) where the supertrace contains an 
evolution operator U of an associated quantum ID prob- 
lem, T is a time-ordering symbol, and . . . stands for op- 
erators that represent the ends of paths and correspond 
physically to density, current, etc. In this form, the av- 
erage can be taken to obtain moments of physical quan- 
tities, and we leave this implicit in later notation. In 
this ID problem vertical rows of links of the network 
become sites, and the vertical direction becomes (imagi- 
nary) time (we assume for the present periodic boundary 
conditions in both directions). The evolution operator 
U , composed of transfer matrices for links and nodes, 
acts in a tensor product of Fock spaces of bosons and 
fermions on each site. The presence of a fermion or bo- 
son on a link — i.e. on a site at an instant of discrete 
time — represents an element of a path traversing that 
link H. Both bosons and fermions are needed to en- 
sure the cancellation of contributions from closed loops. 
Usually one needs two types of bosons and fermions, re- 
tarded and advanced, to be able to obtain two-particle 
properties. However, the particle-hole symmetry relates 
retarded and advanced Green's functions 0. Hence, for 
the study of mean values of simple observables, we need 
only one fermion and one boson per spin direction per 
site. (To study fluctuations and other observables, N 
types of fermion and boson are needed, and the SUSY 
below becomes osp(2JV|2iV) M.) We denote them by / CT , 
b a for the sites related to the links going up (up-sites), 
and f ai b a for the down-sites. On the up-sites, f ai b a are 
canonical, but to ensure the cancellation of closed loops 
we must either take the fermions on the down-sites to 
satisfy {f a , fl,} — —8 aa i, or similarly for the bosons. 

To begin, we consider the spin-rotation invariant case 
with ts-\ = tsi = tg. In this case, for any realization of 
the disorder in the scattering matrices, the transfer ma- 
trices commute with the sum over sites of the eight gener- 
ators (superspin operators) of the superalgebra sl(2|l) = 
osp(2|2), similarly to Ref. f|. The generators for each 
site are constructed as all bilinears in the fermions and 
bosons and their adjoints, which are singlets under the 
random SU(2). These are denoted by [0 B, Q 3 , Q±, 
U±, W±, and have similar expressions for the two types 
of sites. Cancellation of closed loops would only require 
invariance under the gl(l|l) subalgebra generated by B, 
Q 3 , V-, and W+. The larger SUSY that exists when 



ts] = tsi 1S a manifestation of the global SU(2) symme- 
try. 

The transfer matrix describing the evolution on a 
link, after averaging over the random SU(2) matrices, 
projects the states on the corresponding site to a three- 
dimensional subspace of singlets of the random SU(2) m. 
On the up-sites these form the fundamental representa- 
tion 3 of sl(2|l), and we denote them as \m), m = 0, 1, 2. 
Similarly, on the down-sites the three singlet states form 
the representation 3, dual to 3, and we call them |m); m 
is the number of fermions on a site of either type. We find 
that |1) has negative squared norm, (1|1) = —1, while 
the others are positive. Thus, after averaging, we have a 
horizontal chain of sites with alternating dual representa- 
tions on the two sublattices and a discrete-time evolution 
along the vertical direction given by the transfer matrices 
at the nodes, which will be specified below. 

We now consider in detail the node transfer matrix Tg 
on a single node on sublattice S. After the averaging, it 
acts in the tensor product 3(g) 3 for the two sites. Because 
of the sl(2|l) SUSY, we find that it takes the form 



t s = 4^i + (i - 4H ® i- 
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Here the first term contains the projection operator Pi — 
\s)(s\ onto the normalized singlet state \s) = J2 m \m) ® 
|m), while in the second term /, / are the identity opera- 
tors on the two sites (note that I = |0)(0|-|1)(I| + |2)(2|). 
The two terms in Ts represent the two ways to sl(2|l)- 
invariantly couple the in- and out-going states at the 
node, such that the incoming state (in the fundamental 
representation 3) flows out unchanged, turning either to 
the right or the left. They can be represented graphically 
as shown at the top in Fig. pi 

When we multiply the transfer matrices together and 
take the supertrace in the tensor product of all sites to 
calculate the partition function Z — STr U, the result is 
given by the sum of all contributions of closed loops that 
fill the links of the network, weighted by factors of either 
i| or (1 — t s ) for each node. Each loop contributes a 
factor coming from the sum over the three states that 
can propagate around the loop, the supertrace str 1 = 1 
taken in the fundamental 3. It is also clear that Z is 
equal to 1, as it is also before averaging. 

The sum over loops on the links of the network is equiv- 
alent to the bond percolation problem on the square lat- 
tice, as follows. In Fig. g, we shade one-half of the pla- 
quettes of the network in checkerboard fashion. The two 
terms in T$ possible at each node either do or do not 
connect the shaded plaquettes, as indicated by the thick 
undirected line segments. At each A- (respectively, B-) 
node we have a horizontal (vertical) line with probability 
PA = t\ {pb = 1 — t 2 B ). Then on the square lattice formed 
by the shaded plaquettes we have the classical bond per- 
colation problem, and the loops are the boundaries (or 
"hulls") of the percolation clusters. This SUSY repre- 
sentation of percolation easily generalizes to sl(n + l|n) 
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FIG. 2. Graphical representation of the transfer matrices 
Ts, percolating clusters, and two points n, f2, ona hull. 

SUSY, n > 1, using the 2n + 1-dimcnsional fundamental 
representation and its dual. 

Many critical exponents for 2D percolation are known 
exactly. First, there is the correlation length exponent, 
which immediately gives the localization length for the 
spin QH transition, 
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with v = 4/3 [p|; the critical values are pAc = PBc = 1/2 
in the isotropic case. Then, because the basic opera- 
tors of our system are the superspins, which act on the 
states that live on the hulls of the percolation clusters, 
we should consider the exponents associated with these 
hulls. These include an infinite set of scaling dimensions 
for the so-called n-hull operators [M , 



X„ = (An 2 - 1)/12. 
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The exponent x n describes the spatial decay at critical- 



ly ' 



\ri-r 2 \ 



of the probability that n distinct hulls 



each pass close to each of the two points r\ and r 2 ■ There 
is also a set of analogous exponents for the same correla- 
tors near a boundary pT| , 



x n = n(2n — l)/3. 



(5) 



We will now relate further physical quantities within 
our model to percolation exponents, through the SUSY 
mapping. We write the superspins as a single 8- 
component object J for either up- or down-sites. These 
can be inserted in any links of the network, to obtain a 
correlator such as (Q3(ri)Q3(r 2 )), where n and r 2 repre- 
sent links of the network. Then using the same graphical 
expansion, we obtain a sum over loop configurations, now 
with the positions of the insertions marked on the loops, 



and for loops with insertions, the factor 1 is replaced by 
a supertrace (in the fundamental) of the product of ma- 
trices that represent the J's inserted. We then require 
only the total probabilities that loops pass through the 
marked points in various ways. The simplest example is 
a two-point function of J's, which is nonvanishing only 
if the J's are on the same loop, because str J = for all 
components of J. The leading term in the probability 
that the two points are on the same loop (hull) is gov- 
erned by the leading 1-hull operator in the continuum 
theory M], giving 



((-l) il+i2 J(r 1 )J(r 2 ))-|r 1 -r 2 | 
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at the transition, where X\ = 1/4 as specified above. The 
reason for the staggering factors (— l) n , where i\ is the 
site corresponding to r\, will become clear momentarily. 
It is useful to consider the anisotropic limit of the 
model. This is defined by tA, ts —> with a fixed 
ratio tA/ts = (■ Then the transfer matrices Tg 
may be expanded in ts and recombined in the expo- 
nential. The evolution operator has the form U ~ 
exp(— 2tAts JdrTCiu), where the effective Hamiltonian 
7Yid describes a ID superspin chain, with alternating 3 
and 3 representations, and continuous imaginary time r: 
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Here J • J denotes the sl(2|l) invariant product M. The 
transition point, where Hid for an infinitely-long chain 
is gapless, is now at e = 1. The two-site version of Hid 
appeared in Ref. || . 

The sum ^\ Ji is the generator of global SUSY trans- 
formations, and so Ji, viewed as a function of i, is the 
superspin density on the lattice, which gives a subleading 



contribution 



to the J- J correlation at criticality. 



The 1-hull operator must therefore be the staggered part, 
{-lfJi- 

The 1-hull operators represented by (— 1) ! Ji have sev- 
eral physical applications. Components such as Q+ = 
ftf] on the up-sites create fermions, so produce ends 
for the quasiparticle paths. The sum of all such paths 
between n and r 2 represents the quasiparticle Green's 
function, G. To obtain nonzero results on averaging, we 
must multiply the retarded and advanced Green's func- 
tions before averaging, but this can be replaced by a 
spin-singlet combination of our fermions or bosons [pi. 
The staggered part of this averaged correlator represents 
the average zero- frequency density-density ( "diffusion" ) 
propagator |G| 2 (and also the average conductance be- 
tween two point-contacts p3|), which therefore falls as 
|t*x — t^I -1 / 2 at the transition. Moreover, the local den- 
sity of states p(r, E) is represented by another compo- 
nent of the 1-hull operator, because both it and the den- 
sity operator contain wavefunctions squared, ~ l^l 2 ! in 
the original problem and so scale in the same way. The 



energy E itself (set to zero hitherto) has scaling dimen- 
sion yi — 2 — x\ because an imaginary part r\ of E in- 
duces a staggered "magnetic field" term Y^iVPii^v) m 
Tin) . Hence for the average we have at criticality 



p(r,E)~\E\ Xl ' Vl = |£| 1/7 - 



(8) 



Also, since a uniform Zeeman term H causes a shift in 
the Fermi energy, it induces a correlation length £# ~ 
\H\- U1 , where v x = 1/yi = 4/7. 

We have already identified the value v — 4/3 of the 
localization length exponent v with that in percolation. 
In terms of Hid, the effect of a small deviation S = e — 1 
is to add the perturbation <5^i( — l) 4 ^ ■ Ji+i to the 
critical Hid- This term contains the dimer operator 
Di = (— iyji- Ji+i, which is odd under reflection through 
any lattice site (parity). The scaling dimension x 2 of 
the 2-hull operator is the same as that of this "thermal" 
perturbation for the transition, that is v = v 2 = l/$/2, 
2/2 = 2 — X2 |1(| • We therefore expect that the 2-hull op- 
erator is part of a multiplet of staggered two-superspin 
operators, that are similar to Di, but are not all sl(2|l) 
singlets. 

As a final perturbation of the critical Hamiltonian, we 
consider the effect of tsf ^ isj.- This breaks the global 
SU(2) symmetry, and breaks the SUSY to gl(l|l). Taking 
i-Aa = tsa , we find that the effect is to add to Hid a term 
(*T — *l) 2 J2i Ji ' Ji+i> where Jj is the 4-component set of 
generators of gl(l|l), and the product is invariant under 
this algebra. This term is an anisotropy in superspin 
space. The two QH transitions it produces H] cannot 
be seen in our formulation without explicitly introducing 
both retarded and advanced fermions and bosons, and 
we will see only exponentially decaying correlations. The 
correlation length £a induced by A = i-j- — t^ scales as 
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in the notation of Ref. j|], for small A. If the spin 
anisotropy Ji ■ Ji+i has dimension x' , then we will have 
/i = 2/(2 — x'). The operator does not appear to be the 1- 
hull operator, and has the opposite parity to the 2-hull. 
However, the operator product of two 1-hull operators 
has the correct parity and might contain this operator. 
In conformal field theory, the 1-hull operator can be rep- 
resented by 02.2 in the Kac classification of c = Vira- 
soro representations. The fusion rules for this primary 
field with itself contain the leading nontrivial operator 
^1,3, which we view as a subleading 1-hull operator, with 
scaling dimension x\ = 2his = 2/3. We suggest that 
x' = x i = 2/3, which yields \i — 3/2. We further sug- 
gest that this operator describes a random Zeeman term 
(with zero mean). 

Finally, we note that the average two-probe conduc- 
tance of our system with open ends jq], and with tsi = 
tsi, can be related to the number n of hulls that connect 
one end to the other (and back) . Each such configuration 



of loops contributes n to the conductance, times 2 for 
spin, so the mean conductance has the scaling form 



g = 2j2nP(n,L/W,L/0, 



(10) 



n=l 



where P(n, L/W, L/£) is the probability that exactly n 
hulls run from one end to the other and back, for a sys- 
tem of size L by W. This can be considered both for 
periodic and reflecting transverse boundary conditions. 
At the transition, £ = oo, and for large L/W, it is known 
IH that P(n,L/W,0) - e -2^x n L/w for per i dlc, and 
~ e -^x n L/w £ or re fj e cting boundaries. The sum for g is 
dominated by the n = 1 term in this limit, so it has the 



form g 



-L/Zi 



giving the behavior of the localization 



length £id, the only parameter that enters in the com- 
plete distribution of conductance in this limit fll4| . As 
L/W — > 0, we expect that g oc W/L, implying that there 
is a nonzero critical conductivity. 

We may now compare our results with those of recent 
numerical work. In Ref. H, the results obtained were 
v ~ 1.12 and /i ~ 1.45. These are in fair agreement with 
our predictions, especially for \i where our theoretical 
argument is less well established. The authors of Ref. |j| 
study the SUSY spin chain numerically, and find critical 
exponents Xi = 0.26 ± 0.02 and x 2 = 1.24 ± 0.01, in 
excellent agreement with our predictions. 

To conclude, we have used SUSY methods to find a re- 
markable equivalence of a quasiparticle localization prob- 
lem, the spin quantum Hall transition, to 2D percolation, 
resulting in the exact values of three exponents, and the 
universal ratios for the localization length in the ID limit. 
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